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A general class of spreading processes with non-Markovian dynamics 
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Abstract — In this paper we propose a general class of models 
for spreading processes we call the SI*V* model. Unlike many 
works that consider a fixed number of compartmental states, we 
allow an arbitrary number of states on arbitrary graphs with 
heterogeneous parameters for all nodes and edges. As a result, 
this generalizes an extremely large number of models studied 
in the literature including the MSEIV, MSEIR, MSEIS, SEIV, 
SEIR, SEIS, SIV, SIRS, SIR, and SIS models. Furthermore, we 
show how the SI*V* model allows us to model non-Poisson 
spreading processes letting us capture much more complicated 
dynamics than existing works such as Information spreading 
through social networks or the delayed incubation period of 
a disease like Ebola. This is in contrast to the overwhelming 
majority of works in the literature that only consider spreading 
processes that can be captured by a Markov process. After 
developing the stochastic model, we analyze its deterministic 
mean-field approximation and provide conditions for when the 
disease-free equilibrium is stable. Simulations illustrate our 
results. 

1. Introduction 

The study of spreading processes on complex networks has 
recently gained a massive surge of interest. With the wide 
range of applications including the spreading of a computer 
virus, how a product is adopted by a marketplace, or how 
an idea or belief is propagated through a social network, it 
is no surprise that a plethora of different models and studies 
have been devoted to this. However, an overwhelming ma¬ 
jority of the stochastic models proposed and studied assume 
that transitions from one state to another (e.g., a healthy 
individual recovering from a disease) is a Poisson process 
that follows an exponential distribution. Unfortunately, this 
simplifying assumption restricts the applicability of such 
models to exclude a number of specific processes like how 
information is disseminated through Twitter or how the Ebola 
virus is spreading in West Africa. 

In this paper we propose a very general class of epidemic 
models and show how it can be used to study spreading 
processes that don’t necessarily evolve according to an 
exponential distribution. In addition to being able to account 
for non-Poisson spreading processes, our model generalizes 
almost every model studied in the literature including the 
MSEIV, MSEIR, MSEIS, SEIV, SEIR, SEIS, SIV, SIRS, 
SIR, and SIS models. The development and analysis of 
such a general model also allows rapid prototyping of future 
spreading processes that might not even exist today. 


Literature review 

One of the oldest and most commonly studied spread¬ 
ing models is the Susceptible-Infected-Susceptible (SIS) 
model [1]. Early works such as the one above often consider 
simplistic assumptions such as all individuals in a population 
being equally likely to interact with everyone else in the pop¬ 
ulation [2]. One of the first works to consider a continuous¬ 
time SIS model over arbitrary contact graphs using mean 
field theory is [3], which provides conditions on when the 
disease-free state of the system is globally asymptotically 
stable. 

In addition to the simple SIS model, a myriad of different 
models have also been proposed and studied in the literature. 
In [4], [5], the authors add various states to model how 
humans might adapt their behavior when given knowledge 
about the possibility of an emerging epidemic. The work [6] 
considers the possible effect of human behavior changes for 
the three state Susceptible-Alert-Infected-Susceptible (SAIS) 
model. In [7], a four-state generalized Susceptible-Exposed- 
Infected-Vigilant (G-SEIV) model is proposed and studied. 
This model is appealing because it was shown to generalize 
a large number of other models already studied in the 
literature [8], [9]. These models have been used to study 
the propagation of computer viruses [10], [11] or products 
and information [12], and of course the spreading of dis¬ 
eases [13]. However, a large drawback is that all of the works 
above consider an underlying Markov process that drives the 
system. 

While this may be well suited for a number of spreading 
processes, they have also been applied in areas for which 
this is not a very good approximation. A notable example 
is the spreading of the Ebola virus. The work [14] looks at 
the spreading of Ebola in Congo and Uganda in 2004 and 
estimates the spreading properties of the virus fitted to a four- 
state SEIR model. Similarly, the work [15] looks at the more 
recent outbreak of Ebola in West Africa and again estimates 
the parameters of the virus fitted to a six-state model. The 
larger number of states allows the model to better capture 
things like human behavioral changes and also the incuba¬ 
tion period of the Ebola virus. However, just like all the 
works mentioned above, all transitions are assumed to evolve 
according to an exponential distribution. More specifically, 
once an individual is exposed to the virus at some time 
to, the probability that the individual has started showing 
symptoms by time t is given by P{t) = 1 — for 

some £ > 0. However, this is far from a good approximation 
when looking at the empirical data collected over the years. 
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vims and concludes that the incubation period of the disease 
is much better modelled as a log-normal distribution with 
a mean of 12.7 days and standard deviation of 4.31 days, 
which cannot be well captured by the exponential distribution 
above. Another prominent example today is the spreading of 
information through social networks on the internet, such 
as Twitter or Digg. It has been observed multiple times 
that the spreading of information in these networks is again 
better modelled as a log-normal distribution rather than an 
exponential one [17], [18], [19]. 

We are only aware of very few works that have considered 
spreading processes without exponential distributions. The 
work [20] studies the drastic effects that non-exponential 
distributions can have both on the speed of the spreading 
and on the threshold conditions for when the disease will 
die out or persist. A simple SI model is studied in [21] 
without the complexity of an arbitrary graph structure. The 
SIS model with general infection and recovery times is 
considered in [22]. In this work we generalize this idea to a 
much wider class of epidemic models. 

Statement of contributions 


Ebola virus in Section |Vl] Finally, we gather our concluding 
remarks and ideas for future work in Section IVIII 

II. Preliminaries 

We denote by K and IR>o the sets of real and nonnegative 
real numbers, respectively. We define the indicator function 
Iz to be 1 if Z is true, and 0 otherwise. 

Graph theory: Given a directed graph Q with N nodes, 
we denote by A G the associated adjacency matrix. 

The components of A are given by Oji = 1 if and only if 
there exists a directed edge from i to j on the graph Q. 
We denote the in-neighbors and out-neighbors of node i 
as A/]" = {j e W} I = 1} and A/J’"' = {j G 

N} I Oji = 1}, respectively. 

Given a vector q G M", we let diag {qi,..., g„) denote 
the n X n diagonal matrix with gi ,..., on the diagonal. 
Given an arbitrary matrix Q G we define deg(Q) = 

diag >I]"=i the diagonal mxm matrix 

with row sums of Q on the diagonal. For a square matrix Q, 
we define the Faplacian C{Q) = degjQ) — Q. A square 
matrix Q is Metzler if its components qij > 0 for all i f j. 
The following result will be useful in our analysis later. 


The contributions of this paper are threefold. First, we 
propose the SI*V* model that generalizes a very large 
number of models studied in the literature today. Our model 
allows an arbitrary number of ‘infected’ and ‘vigilant’ states 
unlike many models that have a fixed number of states. 
Multiple infected states allows us to capture various stages of 
a disease or spreading process which may have very different 
properties in terms of contagiousness, chance of recovery, 
etc. Multiple vigilant states allows us to capture different rea¬ 
sons that an individual might not be susceptible to a disease 
including behavioral changes, vaccinations, or even death. 
Second, we develop and analyze the deterministic mean- 
field approximation for the model and provide conditions 
for which the disease-free states of our model are globally 
asymptotically stable. Finally, we show how allowing our 
Markov model to have an arbitrary number of states can be 
used to approximate non-Poisson spreading processes with 
arbitrary accuracy. This allows us to much more accurately 
describe real-life phenomena, such as the propagation of 
information through social networks or spreading of the 
Ebola virus, which have recently been shown to evolve 
according to log-normal distributions rather than exponential 
ones. 


Organization 

We begin in Section by reviewing some preliminary 
concepts that will be useful in the remainder of the paper. 
We develop our proposed SI*V* model in Section [In| and 
analyze its stability properties in Section IV In Section 
we show how our general class of models can be used 
to model many stochastic spreading processes with state 
transition times that do not obey an exponential distribution. 
We demonstrate the efficacy of our model and validate our 
stability results through simulations of the spreading of the 


Lemma II.l (Properties of Metzler matrices [23], [24]) 

Given a Metzler matrix Q, the following statements are 
equivalent: 

(i) Q is Hurwitz 

(ii) There exists a positive vector v such that Qv < 0 

(iii) There exists a positive vector w such that nf^Q < 0 

(iv) Q is nonsingular and Q~^ < 0 

III. Model description 

We begin by formulating the SI*V* that we study in 
the remainder of the paper. We follow the idea of the N- 
intertwined SIS model developed in [25] and its extensions to 
the SAIS and generalized SEIV models developed in [6], [7]. 
Although the latter models have been shown to generalize 
many different models studied in the literature, they assume 
a fixed number of compartments for a given disease. Instead, 
we build on a compartmental model studied in [26] in which 
the number of states relating to the disease are arbitrary. 
Unlike the population model (i.e., no graph structure) studied 
in [26], we are interested in proposing and analyzing this 
model applied to arbitrary networks. In Section |V] we show 
how this model can be used to approximate a wide class 
of spreading processes on networks for any type of state 
transitions with arbitrary accuracy, rather than only Poisson 
processes (exponential distributions). 

We consider a virus spreading model with three classes 
of states called the SI*V* model. The first class has only 
one state which is the susceptible state S. The susceptible 
state S corresponds to a healthy individual who is capable 
of being exposed to the disease. The second class has m > 0 
states known as infectious states I. In the infectious class, an 
individual can be in any of the m states given by for k G 
m}. This allows the possibility to model a number of 
variations to the infectious state including human behavior. 



(5*^ Recovery rate from to 
gfcfc Infection internal transition rate from to 
Vigilant internal transition rate from to 
7 I Rate of becoming susceptible from to S 
Rate of becoming vigilant from S to 
0^^ Infection rate due to infected (I^) neighbor j from S to 

TABLE I 

Parameter dehnitions 


severity of the disease, etc. The last class has n > 0 states 
known as vigilant states V. In the vigilant class, an individual 
can be in any of the n states given by for £ G {1 ,..., n}. 
The vigilant class captures individuals who are not infected, 
but also not immediately susceptible to contract the disease. 
The various states in the class can be used to model different 
reasons that the individual is not susceptible such as being 
vaccinated, having just recovered, or even deceased. 

Consider a network with with N nodes. For each node 
i G TV} we define the random variable Xi(t) G 

{S', ..., I"‘, , C"} as the state of node i at a given 

time t. We consider a general directed contact graph Q over 
which the disease can spread. A susceptible node is only able 
to become exposed if it has at least one neighbor that is in 
any of the infectious states. A node i can only be infected 
by nodes in A/}" and can only infect nodes in 

The compartmental Markov process is defined by the 
following parameters. Let be the recovery rate of node 
i going from infectious state to vigilant state V^. This 
allows the possibility to model different recovery rates de¬ 
pending on which state of an infection the individual is in 
and which vigilant state the individual will end up in. We 
let Di = be the matrix that describes these 

transitions. Let be the rate at which an individual in 
infectious state moves to infectious state . This can 
model the various different stages or severity of a disease 
and how individuals move from one stage to another. We let 
Ei = [e’l^'\kk' € be the matrix that describes these 

transitions. We denote by the internal transition rate from 
vigilant state to . This can model different levels and 
types of vigilance in individuals, such as behavioral changes, 
changing medications/vaccines, etc. We let Mi = \ui G 
j^nxn j-jjg matrix that describes these transitions. We 
denote by 7 } and 0 } the rates of moving from to S and 
S' to respectively. Finally, let be the effect that a 
neighbor j G A/]" of node i in state has on i. The rate 
that an individual i moves from S to is then given by 

m 

E E /sAT. 

fc' = l iGTVf 

where 

Y^{t) = lx-(t)=i>‘- ( 1 ) 

Note that when a susceptible individual i becomes infected, 
it always moves into the first infectious state It is then 


Susceptible ' 



Fig. 1. -state SI*V* compartmental model for node i. Internal 

transition rates among the infected and vigilant classes are not shown. 


free to move to the other infectious stages according to Ei. 
All the disease parameters are nonnegative. Table |I] presents 
the definitions of these parameters for convenience. 

The dynamics of the epidemic spread is then modeled 
using the definition of the infinitesimal generator from [27]. 
For brevity, we only write out a small subset of the possible 
transitions, 

m 

P(X,{t') = I^\X,{t)= S,X{t))= ^ ^ Pf.Yf +0, 

P{X,{t') = I^\X,{t) = I\x{t)) = + o, 

P{X,{t') = I^\X,{t) = I\ X{t)) = + o, 

P{XS’) = V^\X,{t) = l\X{t)) = + o, 

where T = f -b At and o = o{At). 

Figure shows the (1 -b m -b n)-state SI*V* compart¬ 
mental Markov model for a single node i. Note that the only 
graph-based transition is from the susceptible state S to the 
first infected state I^. The state of the entire network X{t) 
then lives in a (l-bm-bn)^ dimensional space making it 
very hard to analyze directly. Instead, we utilize a mean- 
field approximation to reduce the complexity of the entire 
system. We do this by replacing in Q by its expected 
value E[Yj]^. 

We denote by [5, (t ), (f),..., I™ (t ), L/ (t),..., L," (f)] ^ 

the probability vector associated with node i being in each 
of these states, i.e., 

m n 

+ = ( 2 ) 
fc'=l i>-l 

for all z S {1,..., N}, k G m}, and £ G nj. 

The SI*V* model we consider in this paper is then given 








by 


m = Y. 7fV/' - d(s. - E E > 

^■'^1 ieA/'f 


’'=1 


int) = s.Y. E /3^^'-^E^“' 

k'-ijeMi" e'=i 


k' = l 


•'=1 


£'^1 k'^1 

n m 

i^t) = -7^ ^ ^ Zf'ef', 

k'^1 

m m n 

vl{t) = E -Y. ^/') 

n 

- + E 


( 4 ) 


'=1 


■'=1 


Next, we are interested in studying the stability properties 
for this set of N{rn + n) ODEs. 


Let X, 


IV. Stability ANALYSIS OF S'/*]/* 

= Jji,..., ir]y. = ., F,"]X = 

[xf,...,x^] , and y = [yf,...,y^] . Naturally, we 
are interested in conditions that will drive the system to 
a disease-free state, i.e., x —)• 0. We begin by writing the 
dynamics of the system as 


X 

y 


= w 


X 

y 


+ H, 


( 5 ) 


where W G ]g^N{n+m)xN{n+7n) Qapl;^J•es the linear part of 
the dynamics and H = with G and 

Hy G M^”, captures the nonlinear part. For convenience, we 
split W into smaller matrices such that 

y^xx y^xti 


w = 


Wyx yVyy 


First, we define the matrix Wxx G that describes 

how infected states affect other infected states. Let 'Wxx = 
[>V*4] be the block matrix where 


W”, = -C{E,) - deg(A) 

describes the internal transitions between the infected states 
and the transitions from all infected states to all vigilant 
states, and 


+ ( 3 ) 

k'^1 

n m 

i^{t) = + E 

k'^i 

m n 

v!{t) = E - ilvl + E ’ 

k' = l i' = l 

for fc G {2,..., to} and 7 G {!,..., n}. 

Due to constraints one of the equations ([^ is redun¬ 
dant. By setting Si{t) = 1 - ^1^=1 (0 “ YTe'=i Vf i't), 

we can describe the system by 

m n m 

il{t) = (1 - E - E ^/') E E 

fe' = l jGTVI" 


Wt = 


« 1 . A™ 

HlJ 5 • • • ? HlJ 

®(m—l)xm 


describes the (linear) transitions of node i from the suscep¬ 
tible state to the first infected state k due to other nodes. 
Second, it is easy to see from 0 that Wxy — OjVm xNn 
because the infected states are not (linearly) affected by the 
vigilant states. 

Third, we define the matrix 'Wyx that describes how 
infected states affect the vigilant states. Since all these 
transitions happen internally (i.e., do no depend on the 
network structure), Wyx = \Wyx\ G ^ block 

diagonal matrix where 


W” = - 

^ ^ yx % 


d} 

Of 

6 »r 


Of 

el 


describes the transitions from all infected states to all vigilant 
states and the transitions from the susceptible state to all 
vigilant states. 


Finally, we define the matrix >Vj,y G 


'JnxNn 


that 

describes how vigilant states affect other vigilant states. As 
before, since all these transitions happen internally, Wyy = 
[W”y] is a block diagonal matrix where 


= -C{Mi) - 


Of 

Of 

e? 


Of 

Of 

0 " 


-diag(7},...,7'") 


describes the internal transitions between the vigilant states, 
transitions from the susceptible state to all vigilant states, and 
transitions from all vigilant states to the susceptible state. 

We now define the column vector H. As can be seen in 
equation 0, the nonlinearities only enter into the dynamics 
of the first infectious state k. Thus, we can describe HI G 


/ Tk' _ Y^ 

Z^fc' = l = l ''i J Z^/c' = l Z^, 


■jeHi 


Hij 


Of" rr),— 


(m—l)xl 


Due to removing the susceptible state, we also have a 
constant forcing given by H^ = [9f,..., 0”] ^ G K". 

Theorem IV.l (Sufficient condition for global stability of 
disease-free equilibrium) The disease-fr^^ states of 0 are 
globally asymptotically stable if Wxx Is Hurwitz. 
















Proof: We begin by noticing that Wxx is a Metzler 
matrix. Using Lemma II. 1 and the condition that Wxx is 
Hurwitz, we know there exists a positive vector V G 
such that v^Wxx < 0. Consider the Lyapunov function 

J = 


then 


j = v'^x 

= v'^Wxx^ + v'^Hx 
< v'^Wxx^ 


is strictly negative for all x ^ 0. It can then easily be shown 
using LaSalle’s Invariance Principle that x —> 0 [28]. ■ 

Note that this conservative result essentially ignores the 
vigilant class and how it can help a disease die out. In other 
words, this result states that if the vigilant class is removed 
and the disease-free state is still globally asymptotically 
stable, adding the vigilant class cannot hurt. For this reason, 
we do not see any forms of the parameters , 7 ^, or 9f 
appear in the condition. 

In order to find the necessary and sufficient condition, 
we linearize the entire system around the disease-free equi¬ 
librium. For simplicity, we assume here that there exists 
no absorbing state inside the vigilant class. Note that if 
there exist any absorbing states in the vigilant class, it does 
not make sense to discuss asymptotic properties anyway, as 
all individuals will eventually end up there, resulting in a 
disease-free equilibrium. 

We begin by computing the unique equilibrium for the 
vigilance class states. Letting y = 0 and x = 0, we get 

y* = -Wfy^Hy. 

The inverse of Wyy is guaranteed to exist because it is a 
block diagonal matrix made up of negative definite matrices. 
Let [U/, = y* for all i G N}. The 

linearization of ([5]l around the point x = 0 , y = y* is then 
given by 


X 

= Q 

X 

_ y 


y 


where 

- deg(A), 


= 1 - E 


"=1 


Pij 


■'ij 

0 


(m—1) X m 


Qxy — yVxy — 0 


NmxNn 


5 Qyx - '^^yxi Qyy - yVy 


Theorem IV.2 (Local stability of disease-free equilib¬ 
rium) The disease-free states of 01 are locally asymptot¬ 
ically stable if and only if Qxx is Hurwitz. 

Proof: Since the linear dynamics of the infectious 
states X do not depend on the vigilant states y, we can write 

^ — Qxx^. 



Exposed I 



Infected i 



I 

I 

I 


Fig. 2. Four-state G-SEIV compartmental model for node i. 


It is then clear that x —^ 0 in this linearization of the original 
system if and only if Qxx is Hurwitz. ■ 

Remark IV.3 (Global stability of disease-free equilib¬ 
rium) We conjecture, and verify through simulation, that the 
necessary and sufficient local stability result of Theorem |IV.2| 
is indeed a global result but it has yet to be shown. We expect 
to have its proof completed in the near future. • 

We note here that based on these results, determining the 
stability properties of the disease-free equilibrium amounts 
to checking if an Nmx Nm matrix is Hurwitz, even though 
the original system is N{m n)-dimensional. 

V. Non-exponential distributions 

In this section we show how existing epidemic models 
can be studied with non-Markovian state transitions by 
appropriately constructing an instance of our SI*V* model. 
To simplify the exposition, we demonstrate this idea for a 
single example, the spreading of the Ebola virus, but note that 
the same idea can be used to expand a very large number of 
different models. 

The underlying model we consider is the four-state G- 
SEIV model proposed in [7] and shown in Eigure The 
‘Susceptible’ state S captures individuals who are able to 
be exposed to the disease, the ‘Exposed’ state E captures 
individuals who have been exposed to the disease but have 
not yet developed symptoms, the ‘Infected’ state I captures 
individuals who are displaying symptoms and contagious, 
and the ‘Vigilant’ state V captures individuals who are not 
immediately susceptible to the disease (e.g., just recovered, 
quarantining themselves, deceased). 

We assume that all transitions are Poisson processes except 
the transition from exposed to infected. Eor this transition 
we use a log-normal distribution with mean 12.7 days and 
standard deviation 4.31 days as proposed in [16] based on 
empirical data. In the constraction of the approximation, 
an important role is played by the class of probability 
distributions called phase-type distributions [29]. 

















Consider a time-homogeneous Markov process in 
continuous-time with p -f 1 (p > 1) states such that the 
states 1, ..., p are transient and the state p -|- 1 is absorbing. 
The infinitesimal generator of the process is then necessarily 
of the form 

S -51„xi' 

/ , (6) 

Oixp 0 

where S G is an invertible Metzler matrix with non¬ 

positive row-sums. Also let 



0 5 10 15 20 25 30 35 


(7) 

be the initial distribution of the Markov process. Then, the 
time to absorption into the state p-|-1 of the Markov process, 
denoted by {4>,S), is called a phase-type distribution. It is 
known that the set of phase-type distributions is dense in the 
set of positive valued distributions [30]. Moreover, there are 
efficient fitting algorithms to approximate a given arbitrary 
distribution by a phase-type distribution [29]. 

We now show how this can be used to expand the G- 
SEIV model to an instance of our SI*V* model such that 
the time it takes from to reach the infected state from the 
exposed state follows a phase-type distribution. 


</> 

0 


Proposition V.l Consider the SI*V* model with m = p-fl 
infectious states, where corresponds to the infected state 
and for k G {1 ,... ,p} correspond to the exposed state. 
Let Si G be given by 


[Sl]kk' 


k f k', 

— Yl'k'=i Otherwise. 


( 8 ) 


Then the length of time that it takes a node i to go from 
state to follows the distribution (ei,S'/), where ei = 

[l,0ixp]^. 


Proof: Without loss of generality, assume that a node i 
becomes exposed at time t — Q, i.e., 2fi(0) = I^. Letting T 
be the time that node i first reaches the infected state, i.e., 
Xi(T) = I™, we need to show that T follows (ei,5/). 

We begin by noticing that the initial condition -Ai(0) = Ci 
corresponds to node i having just entered the exposed state. 
This is because in our SI*V* model, a susceptible node 
can only enter the infected states through I^. Also, from the 
definition of the model, it is clear that the Markov process Xi 
has the infinitesimal generator on the time interval [0, Tj. 
Also /"* is clearly the only absorbing state of the process 
(when confined on [0, Tj). The above observation shows that 
T follows the phase-type distribution (ei,5/) because we 
have satisfied the definitions of the infinitesimal generator (|^ 
and initial distribution 0- ■ 

Proposition |V.1| shows that it is possible to model the 
transition from the exposed state to the infected state of 
the SEIV model as a phase-type distribution. This is done 
by essentially expanding the exposed state in the original 
SEIV model from a single state to p. As noted earlier, this 


Fig. 3. Approximation of log-normal distribution (mean 12.7 days and 
standard deviation 4.31 days) with phase-type distributions for Ebola’s 
incubation period [16]. 


is only one specific example that can be extended to model 
many different state transitions as phase-type distributions 
rather than exponential ones. Next, we show how an arbitrary 
distribution can be approximated as a phase-type distribution 
and how to choose the parameters for our SI*V* model to 
realize the desired distribution. 

Continuing with our Ebola example, we show how phase- 
type distributions can well approximate the log-normal dis¬ 
tribution of Ebola’s incubation time with mean of 12.7 days 
and a standard-deviation of 4.31 days. To do this, we utilize 
the expectation-maximization algorithm proposed in [29] 
Eigure shows the results for different amounts of internal 
states p. We can see here that the more internal states p we 
use, the closer our phase-type distribution becomes to the 
desired log-normal distribution. An instance of the phase- 
type distribution for p = 10 is shown in Figure 


Remark V.2 In fact one can show that the set of the phase- 
type distributions of the form (ei,iS) is dense in the set of 
all the phase-type distributions as follows. Let {4>,S) be a 
given phase-type distribution. For a positive real number r 
define the block-matrix 


Tr 


—r rcjf^ 
0 5 


(9) 


Then we can prove that the sequence of the phase type 
distributions {(ei, with p2 states converges to 

the phase type distribution (</>, 5) with p -\- 1 states. Here 
we provide only an outline of the proof. Let X be the 
time-homogeneous Markov process having the infinitesimal 
generator 0 and the initial distribution Ci. Define ti = 
supjr : X{t) = 1}, i.e., the time at which the Markov 
process leaves the first state, and let t 2 be the time the process 
X is absorbed into the state p -\-2. By the definition of the 
first row of the generator T^, we see that X{ti) follows the 
distribution [O, O]^ G 

Therefore <2 —G follows (0,5). Moreover, as r increases, 
the probability density function of G converges to the Dirac 
delta on the point 0. Therefore the random variable ^2 = 
— which follows (ei,Tr) by definition, converges 
to {(j),S). The details are omitted due to space restrictions.* 













The implications of Remark V.2 are that although our 
SI*V* model only allows a susceptible node to enter the 
infected class through I^, we are still able to model any 
phase-type distribution rather than just (ei,5). 


VI. Simulations 

Here we demonstrate how the results of Section |V] can be 
used to model a spreading process with a non-exponential 
state transition and validate the stability results of Section [TV| 
by simulating the spreading of Ebola. The underlying model 
we use is a four-state G-SEIV model proposed in [7] 
and shown in Eigure The ‘Susceptible’ state S captures 
individuals who are able to be exposed to the disease, 
the ‘Exposed’ state E captures individuals who have been 
exposed to the disease but have not yet developed symptoms, 
the ‘Infected’ state / captures individuals who are display¬ 
ing symptoms and contagious, and the ‘Vigilant’ state V 
captures individuals who are not immediately susceptible to 
the disease (e.g., just recovered, quarantining themselves). 

We assume that all transitions are Poisson processes except 
the transition from E to L Eor this transition we use a 
log-normal distribution with mean 12.7 days and standard 
deviation 4.31 days as proposed in [16] based on empirical 
data. Using the EM algorithm proposed in [29] with p = 10 
phases, we expand the exposed state from a single state to 
10 states. This means we can describe our four state non- 
Poisson SEIV model by a 13-state Poisson SI*V* model 
with one susceptible state, one vigilant state, and m = 11 
infectious states where is the only contagious state. The 
other infectious states for k G {1 ,..., 10} correspond to 
the exposed (but not symptomatic) state of the original SEIV 
model. The obtained internal compartmental model between 
the exposed state and infected state is shown in Eigure 

Eor our simulations we consider a strongly connected 
Erdos-Renyi graph A with = 20 nodes and connection 
probability 0.15. Initially, the vaccination rates 6i are ran¬ 
domly chosen from a uniform distribution 9i G [0.3,0.8] and 
the rates of becoming susceptible ji G [0.2, 0.7]. Since it is 
known that Ebola can only be transmitted by people who are 
showing symptoms, we set = 0 for all fc S {1 ,..., 10}. 
Eor links that exist in the graph A we randomly set the 
infection rate fiA G [0.1,0.4]. Similarly, we assume that 
only infected individuals can recover, thus we set = 0 
for all k G {1,..., 10} and randomly set the recovery rate 
Si G [0.1,0.4]. Since we only have one vigilant state, there 
are no internal transition parameters p. 

To demonstrate the effectiveness of the expansion of our 
model to capture the log-normal incubation times of Ebola, 
we randomly set the initial conditions of being exposed to 
li(0) G [0.25,0.75] and the susceptible state to Si(0) = 
1 — 7/(0). Thus, we assume that there are initially no nodes in 
the vigilant V or infected states for k G {2,. ..,11}. Eor 
the parameters used, we obtain Xmax{Qxx) = —0.1264 as the 
largest real part of the eigenvalues of Qxx- Figure a) shows 
the evolution of the maximum, minimum, and average prob¬ 
abilities of being in any of the 11 infected states over time 
F’i(f) = Figure l^b) shows the probabilities 


of being in only the truly infected (and symptomatic) state 
for all nodes i. Here we can see the effectiveness of 
our expanded model in capturing the log-normal incubation 
times of Ebola, seeing the peak of infections at 12.7 days. 
Given enough time, we see that all infections eventually die 
out as the stability condition of Theorem IV.2 is satisfied. 



Days 

(a) 



Days 

(b) 


Fig. 5. Plots of (a) the minimum, maximum, and average trajectories of the 
probability of each node being in any infected state and (b) the trajectories 
of the probability of each node being in the infectious state 

In Eigure we vary the recovery rates Si and infection 
rates (SA and look at the steady-state probabilities P{oo) 
of each node being in any infectious state where we ap¬ 
proximate P{oo) « P{T) for large T. We can see here the 
sharp threshold property that occurs as Amax(Qa;x) moves 
from negative to positive, validating our stability results of 
Section |lV] 


VII. Conclusions 

In this work we have proposed a general class of stochastic 
epidemic models called SI*V* on arbitrary graphs, with 
heterogeneous node and edge parameters, and an arbitrary 
number of states. We have then provided conditions for when 
the disease-free equilibrium of its mean-field approximation 
is stable. Eurthermore, we have shown how this general class 
of models can be used to handle state transitions that don’t 
follow an exponential distribution, unlike the overwhelming 
majority of works in the literature. We demonstrate this 
modeling capability by simulating the spreading of the Ebola 
virus, which is known to have a non-exponentially distributed 











Fig. 4. Approximation of Ebola incubation time distribution with phase-type (p = 10) distribution. 



Fig. 6. Plot of the minimum, maximum, and average steady-state 
probabilities of each node being in any infected state. 

incubation time (i.e., time it takes to show symptoms once 
an individual is exposed). For future work we are interested 
in studying how to control this model which can be used in 
a much wider range of applications than before due to its 
capabilities in modeling non-Poisson spreading processes. 
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